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We show for the first time that the quenching of electronic excitation from nuclear recoils in 
liquid xenon is well-described by Lindhard theory, if the nuclear recoil energy is reconstructed using 
the combined (scintillation and ionization) energy scale proposed by Shutt et al. We argue for 
the adoption of this perspective in favor of the existing preference for reconstructing nuclear recoil 
energy solely from primary scintillation. We show that signal partitioning into scintillation and 
ionization is well-described by the Thomas-Imel box model. We discuss the implications for liquid 
xenon detectors aimed at the direct detection of dark matter. 
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I. INTRODUCTION 

There is considerable experimental effort dedicated to 
the direct detection of particle dark matter. Among the 
various detection strategies PQ, dual-phase liquid xenon 
detectors have recently achieved a rapid scale-up in tar- 
get mass [2J, which directly improves the sensitivity of 
the dark matter search. The expected experimental sig- 
nature from the scattering of a dark matter particle is 
in most cases a low-energy O(keV) recoiling nucleus. A 
fraction /„ of the nuclear recoil energy E„ r is transferred 
to measurable electronic excitation, and the rest is lost to 
heat, e.g. atomic motion. The theoretical prediction for 
the quenching /„ obtained by Lindhard et al. [3] agrees 
to within ±20% with the measured amount of ioniza- 
tion in germanium detectors [H [5J . However, it has been 
observed that the Lindhard theory does not agree with 
measurements of ionization quenching of nuclear recoils 
in liquid xenon [BJ. 

Experiments with liquid xenon as the target medium 
historically use only the primary scintillation signal to 
reconstruct E nr [2] I7HT0] . Commensurate with this non- 
linear energy scale determined by primary scintillation, 
there has been a substantial experimental effort aimed 
at measuring (see [HI [12] and references therein) and 
understanding [13HI6] the scintillation quenching £ e // 
of liquid xenon for nuclear recoils. In spite of these ef- 
forts, factor of two disagreement persists between recent 
measurements, and Lindhard theory does not correctly 
predict £ e // [SI EH (HI EO HE]- However, with semi- 
empirical modifications, reasonable agreement may be 
obtained [TT1 HH HI] - 

Meanwhile, it has been known for nearly a decade [19] 
that a linear, drift-field independent energy scale with a 
substantial improvement in energy resolution [201 121) is 
obtained from a simultaneous measurement of the num- 
ber of primary scintillation photons n 7 and ionized elec- 
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trons n e . This is given for electron and nuclear recoils by 
Shutt et al. [22] as 

E er = e(n 7 + n e ), (1) 

E„r = e(n 7 + n e )/f n , (2) 

where e = 13.8±0.9eV [H] is the average energy re- 
quired to create a single quanta, either photon or elec- 
tron. Other measurements have found e = 13.7 ± 0.2 eV 
[2"3"] and e = 14 eV [24- . Eq. [T] is widely used to measure 
E er for electron recoils from gamma and beta particles 
[HI EDI 123], but Eq. [2] as a measure of E„ r has been largely 
ignored. We note that this formulation is analogous to 
the ionization energy scale in germanium detectors. 

In this article, we show that the Lindhard theory is 
consistent with the quenching of the total (n e + n 7 ) elec- 
tronic excitation from nuclear recoils in liquid xenon. 
This further motivates the adoption of the combined en- 
ergy scale (Eq. [2]), but it does not resolve the disagree- 
ment between recent measurements of C e f / [TTJ [12] . We 
therefore pinpoint the likely source of the present dis- 
agreement, and suggest a strategy for reducing system- 
atic effects in future measurements of either C e ff or /„. 

Our motivation for this work is to obtain the best 
reconstruction of the detected nuclear recoil energy in 
liquid xenon dark matter search experiments. The pre- 
dicted nuclear recoil energy spectrum of halo-bound par- 
ticle dark matter elastically scattering from a xenon tar- 
get falls sharply with increasing E„ r [55]. This implies 
that the sensitivity of a given exposure (detector mass 
x time) of a liquid xenon detector has a strong depen- 
dence on the detector energy scale. For light (< 10 GeV) 
dark matter particles, the dependence on energy scal- 
ing and threshold is even more severe. Recent results 
from XENON100 [2J, in which no candidate events were 
observed, gave a clear example of this. Depending on 
assumptions about the energy scaling (obtained therein 
from primary scintillation via C e ff) and the energy res- 
olution, XENON100 may or may not exclude parameter 
space consistent with potential signal in the CoGeNT [26 
and DAMA [27] experiments. 
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II. RECONSTRUCTING NUCLEAR RECOIL 
ENERGY 

Lindhard et al. [3J calculated a general expression 
for the expected fraction of nuclear recoil energy that 
is transferred to electrons. It can be written as 

f n = k-g/(l + k-g), (3) 

with k = 0.133 Z 2 / 3 ' A~ x l' 2 . Physically, k is a propor- 
tionality constant between the electronic stopping power 
dE/dx and the velocity of the projectile (which in this 
context is a recoiling xenon atom). The relation is most 
simply expressed in terms of dimensionless variables, 
as in [3]. For xenon, Lindhard's calculation results in 
k = 0.166. Recently, Hitachi calculated from first prin- 
ciples the electronic stopping power of recoiling xenon 
atoms in a liquid xenon target. The result is shown in 
Fig. 5 of [6], and discussed further in [16]. In terms of 
the dimensionless variables of [3], his calculation corre- 
sponds to k = 0.110. Note that no analytic form was 
given for the energy-dependent function g in [3J, and we 
have used the parameterization given in [25]. In Fig. [I] 
we show /„ as calculated from Eq. [3] for these two val- 
ues of k (solid and dashed curves). In [3J, Lindhard et 
al. cautioned that "Maybe the greatest uncertainty is 
the proportionality factor, k... [which] is often on the 
interval 0.10 < k < 0.20." Ideally, the remedy for this 
uncertainty may be obtained from data. 

In order to compare with data, we write Eq. [2] as 

SI S9 

U = e(— + — )/E„ r) (4) 

with n 7 and n e in terms of the experimentally mea- 
sured quantities SI and S2. These are just the num- 
ber of recorded photoelectrons in the primary scintilla- 
tion and ionization (measured from proportional scintil- 
lation) signals in a dual-phase xenon detector. The num- 
ber of primary scintillation photons is rt 7 = Sl/a\, where 
ai ~ 0(0.1) is the total efficiency to convert a scintilla- 
tion photon to a detectable photoelectron. The number 
of ionized electrons is n e — S2/a2, where ct2 ~ O(10) 
is the number of photoelectrons registered from the pro- 
portional scintillation resulting from a single ionized elec- 
tron. While o?2 is reasonably easy to measure in dual- 
phase liquid xenon detectors, a± is difficult to measure 
directly. As a result, experiments instead quote the scin- 
tillation light yield L y (units of photoelectrons / keV) 
of a mono-energetic gamma source. The proportionality 
constant between ot\ and L y depends strongly on both 
the incident gamma energy and the electric field (Ed) 
applied across the liquid xenon target. For a 122 keV 
gamma from 57 Co at Ed — 0, a detector-independent 
expression for the proportionality is oi\ — 0.015L y |23j . 

We need a.\ in order to use Eq. |4j and first cross-check 
the relationship given above, which implies ot\ = 0.080 
[23J for the XENON10 detector [SJ. XENON10 mea- 
sured ai — 24 ± 1 [55] from the background distribution 



of single electrons. The value of a\ is also uniquely de- 
termined by requiring Eq. [I] to reproduce the correct 
peak positions of gamma lines. Since the energy scale 
is linear, any mono-energetic source will suffice. From 
the 164 keV gamma observed by XENON10, we find 
ai = 0.078 ±0.005. 
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FIG. 1. Quenching of electronic excitation from nuclear re- 
coils in liquid xenon: from [11] (#), from [11] as corrected 
by [29] uncertainty not shown but similar to •), from 
[12] (O) and from [28] (()). Also shown are the theoretical 
prediction [3] for two calculated values of k (solid and dashed 
curves). 

Simultaneous measurements of the scintillation (C e ff) 
and ionization (Q y ) yield of liquid xenon as a function of 
E nr were obtained by Manzur et al. [11] . Their ionization 
yield data was presented in terms of n e (so we do not need 
02), and their scintillation data in terms of L y for 122 keV 
gammas. Using the scaling relation given above, we infer 
u\ for their detector. We then use Eq. [4] to cast the 
results from [11] in terms of /„. This is shown in Fig. [I] 
(as #. with lcr uncertainty). The combined £ e // and Q y 
measurements of Manzur et al. are not quite consistent 
with the XENON10 nuclear recoil band measurement [7J. 
In [29] it is argued that for the three data points below 
E nr ~ 6 keV, the most likely origin of the disagreement 
is that Q y was overstated by about la due to spurious 
threshold effects. The Manzur et al. data as corrected 
by [29] is also shown in Fig. [l](as -k, uncertainty similar 
but omitted for clarity). 

The experiments described in pi] and [TT] (and, for 
that matter, [5]) obtained their data with different val- 
ues of Ed- Although the values of Ed ranged from about 
0.5 kV/cm 2 to 1.0 kV/cm [11], the effects of this dif- 
ference are negligible. As shown in Fig. 3 of [B], the 
scintillation and ionization signals from nuclear recoils 
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exhibit no significant dependence on Ed, in the range 
0.5 - 1.0 kV/cm. At values of E d < 0.2 kV/cm, a de- 
crease of about 8% is observed in the ionization signal, 
relative to Ed > 0.5 kV/cm. In the same range, the 
scintillation signal is observed to increase by a similar 
amount. 

Aprile et al. recently measured £ e // &t E d = 
[12) . This condition prevented them from making a si- 
multaneous measurement of the ionization yield, as in 
[IT]. Following [29], we used the constraint from the 
XENON10 nuclear recoil band to infer the ionization 
yield (at Ed — 0.73 kV/cm) for this data, taking care 
to account for the 8% shift in n e and n 7 . The resulting 
/„ is shown in Fig. [I] (as O, with la uncertainty). Un- 
certainty in the ionization yield was assumed to be 20%, 
which is similar to the uncertainty reported in The 
slightly smaller total uncertainty in this case results from 
the higher ot\ obtained in that experiment. 

A final data point in Fig. [I] (0) was obtained from 
the endpoint of the XENON10 neutron calibration [25] . 
The maximum recoil energy imparted by Am-Be neu- 
trons to a xenon target via elastic scattering is E nr = 
4E„m n m/(ro n + to) 2 . The maximum incident neutron 
energy was E„ = 10.5 ± 0.5 MeV [3D], with the uncer- 
tainty represented by the horizontal error bar. 



III. DISCUSSION 

Signal quenching from nuclear recoils in liquid xenon 
appears well-described by the Lindhard model of quench- 
ing, once account is taken of all energy transferred to elec- 
trons. This is a very significant result, but it should not 
be surprising. Lindhard et al. defined fj as "the sum total 
of the energy given to electrons," [3J with explicit men- 
tion of both excitation and ionization. In liquid xenon, 
the former leads to n 7 and the latter to both n e and n 7 
[TP] , In our notation, which follows that of other recent 
work, /„ = fj/E nr . 

Considering the uncertainties in the existing data, it is 
unclear which calculated value of k is preferred. But it is 
clear that systematic uncertainty in future measurements 
of f n can be reduced. Essential steps not taken in pre- 
vious work include direct in-situ measurement of ot\ and 
a 2 . A particularly apt choice for measuring ct\ would be 
the internal, homogenous 40 keV line from 83m Kr [3Tll32j . 
Non-uniformities in signal collection at the detector edges 
could be minimized by x — y position reconstruction, and 
rejection of edge events. As pointed out in [2pJ, it is also 
desirable to obtain a measurement of the nuclear recoil 
band, as a cross-check. However, this may be difficult to 
implement without shielding from external background 
radiation. 

It may be preferable to use Eq. [2] to reconstruct nu- 
clear recoil energy, instead of E„ r = SI/ (L y C e f f)S e / S n 
[21 LZl [9] . In the latter, S e and S n are the scintillation re- 
duction due to Ed, for electron and nuclear recoils. In 
either case, it will be essential to retain a distinct mea- 



sure of the primary scintillation response. This is be- 
cause ai <C ct2, and so the low-energy signal detection 
efficiency and energy resolution depend almost entirely 
on SI. Additionally, some experiments rely solely on the 
SI signal |33j . The primary scintillation response can be 
related to /„ using the photon fraction n 7 /(n 7 + n e ), 
as suggested in [33]. In Fig. [2] we show the photon 
fraction (large symbols) and the ionized electron frac- 
tion n e /(n 7 + n e ) (small symbols) for the data from Fig. 
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FIG. 2. The photon fraction n 7 /(n 7 + n e ) (large symbols) 
and the electron fraction n e /(n 7 + n e ) (small symbols) for 
nuclear recoils in liquid xenon, for Ed = 1.0 kV/cm. Symbols 
correspond to data shown in Fig. [TJ Uncertainty in O and ~k 
were omitted for clarity, and are similar to uncertainty in #. 
Curves are a best fit of Eq. [6] as summarized in Table [I] 



It is well-understood that energy loss in liquid noble 
gases leads to a number N ex of excited atoms and a num- 
ber Ni of ionized atoms. We make the usual assumption 
that each excited xenon atom leads to the creation of one 
scintillation photon, and that each ionized atom leads to 
a single electron unless it recombines |19) . Ionized atoms 
which rccombinc result in a single scintillation photon, 
giving Ni + N ex = n e + rt 7 independent of recombina- 
tion. The fraction of electrons collected (i.e. those es- 
caping recombination) is predicted by the Thomas-Imel 
box model [35] to be 



n e 1 . . Na 
Ar£ ln(1 + a ^4oV 



(5) 



which describes Ni initial electron-ion pairs in a box of 
dimension a. This model has been shown to work well 
for spatially small tracks (nuclear recoils and low-energy 
electron recoils) when the term a/ (a 2 v) is held constant. 
[2"3] . The electron fraction F e — n e / (n e + n 7 ) is then 
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given by 



£ n e + n 1 



(6) 



and the photon fraction is 1 — Considering Eq. [2j we 
also have 



(7) 



We fit Eq. |6]to the data, treating a/(a 2 v) and N ex /Ni as 
free parameters and using y* n given by Lindhard with k — 
0.166. The results are shown in Fig. [2] and summarized 
in Table U 



TABLE I. Result of fitting Eq. [6] to the data shown in Fig. 
[2] The first (second) row corresponds to the solid (dashed) 
curves in Fig. [2] The XENON10 data point (0) is not in- 
cluded in the fits. 



data 


N ex /Ni 


a/(a 2 v) 


Ref. 


• 


0.86 


0.028 


ED 




1.05 


0.025 


CD 


★ 


1.04 


0.030 


EDES] 


O 


1.13 


0.042 


EU 



a E„ r > 7 keV only 

We find that for nuclear recoils in liquid xenon, 
N ex /Ni ~ 1 with an uncertainty of about 15%. This is in 
good agreement with [23], which found N ex /Ni = 0.89. It 
has been known for some time that for electron recoils, 
N ex /Ni = 0.06 [!)]. As suggested in [53], the roughly 
factor xl5 difference in initial exciton to ion ratios pro- 
duced by electron and nuclear recoils appears to be the 
origin of the distinct S2/S1 bands for charged versus neu- 
tral particle interactions. We note that Eq. [6] provides a 
less than satisfactory fit to data (•) with E„ r < 7 keV, 
and so we also show the results of a fit which excludes 
this region. The disagreement may indicate that N ex /Ni 
is energy-dependent in this regime. A detailed discus- 
sion of the recombination physics contained in the ratio 
a/(a 2 v) is beyond the scope of this work, and we simply 
note that our results are very similar to those obtained 
in p23J, which used the S2/S1 band centroid rather than 
fixed-energy data points. 

The observed trend in the electron fraction data (• 
and implies that by about E„ r ~ 2 keV, the measured 
signal may consist almost entirely of electrons. If the 
photon fraction is 1 — T e = 0.2 at 2 keV, and we conser- 
vatively assume that Hitachi's calculation of k = 0.110 is 
correct, the expected SI signal at this energy is ~ 3 pho- 
tons. Considering typical values of a\, the probability 
that this results in a measurable number of photoelec- 
trons is nearly zero. From the relation 



£*ff ~~ (1 — -^e) 



fn gl S e 



(8) 



in which all energy dependence is encoded in J- e and /„, 
we would expect £ e // = 0.025. Meanwhile, an electron 



fraction F e — 0.8 implies an S2 signal of ~ 14 elec- 
trons at 2 keV. This translates to an ionization yield 
Q y = 6.8 electrons/keV, or about higher than the 
value obtained in |37| . On the other hand, if the trend 
predicted by Eq. [6] is correct and 1 — T e = 0.5 at 2 keV, 
we would expect £ e ff — 0.063. This implies an S2 signal 
of 8.5 electrons, and the ionization yield Q y = 4.3 elec- 
trons/keV is about ler below the value obtained in [37]. 

It has been pointed out that in order for electronic exci- 
tation to result from a two-body collision between a pro- 
jectile (the recoiling atom) and an atomic electron, the 
maximum possible energy transfer in the collision must 
exceed the band gap energy E ff of the target [38] . This 
condition can be expressed in terms of the projectile ve- 
locity v and the electron velocity v e as E g < 2m e v(v + v e ) 
|40) . The electron velocity is often identified with the 
Fermi velocity [38] [39] v F = (3?r 2 / o e ) 1 / 3 fi/m e . In the ap- 
plication of the Fermi-gas model to metals, the electron 
density p e is calculated only for the valence electrons [H] . 
This follows from the assumption of a "free and indepen- 
dent electron gas." For insulators or semiconductors, the 
situation is less clear. If all 54 electrons in a xenon atom 
are considered, the requisite projectile velocity which sat- 
isfies the inequality (with E g = 9.3 cV in xenon [43]) im- 
plies a cutoff in electronic excitation below a nuclear re- 
coil energy E„ r = 39 keV [H] . We suggest that it is more 
physical to calculate vp under the assumption that the 
Fermi level Ep lies at the midpoint of the energy gap [42] . 
For xenon this implies an effective 3 valence electrons per 
atom, with a predicted kinematic cutoff in electronic ex- 
citation below E nr = 157 keV. Considering Fig. [l] we 
do not see significant evidence for any kinematic cutoff 
in excess of the Lindhard prediction. Applying similar 
reasoning to germanium, one would consider an effective 
1/40 valence electrons per atom in calculating vp, and 
predict a kinematic cutoff below E„ r = 6.5 keV. Again, 
no evidence for such a cutoff is observed 5 . A simple 
interpretation of this is that direct two-body collisions 
are not the dominant electronic energy-loss mechanism 
for recoiling atoms. Lindhard et al. indirectly suggest as 
much [3J. 



IV. SUMMARY 

We have shown that the quenching of electronic exci- 
tation from nuclear recoils in liquid xenon appears well- 
described by the prediction of Lindhard et al. [3J when 
taking the appropriate measure of electronic excitation, 
i.e. the combined energy scale given by Eq. [2] Further 
measurements with reduced uncertainty will be needed 
to elucidate the correct value of k in Eq. [3] and we have 
provided three specific suggestions for achieving this. We 
have urged the adoption of the combined energy scale for 
nuclear recoils, with the photon and electron fractions de- 
scribing the signal partitioning. This has the distinct ad- 
vantage of allowing comparison with the quenching the- 
ory of Lindhard et al, and the recombination model of 
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Thomas and Imel. 

For direct detection experiments with null results, the 
choice of energy scale has a negligible effect on exclu- 
sion limits. However, one eventually expects to observe 
nuclear recoils from the scattering of dark matter — op- 
timistically, in the 100+ live days of blinded data so 
far accumulated by the XENON100 collaboration. At 
that juncture it will be essential to obtain the most ac- 
curate reconstruction of the true recoil energy. We il- 
lustrate this point with a simple example, considering 
the recent XENON100 result (Fig. 3 of 0). Sup- 
pose there were two hypothetical nuclear recoil events 
A and B, both with SI = 4 photoelectrons, but with 
event A having logio(S2/Sl) = 2.3, and event B having 
logi (S2/Sl) = 1.9. The Sl-only energy scale employed 
in [5] assigns the same E nr = 8.7±4.3 keV to both events. 
In contrast, the combined energy scale (Eq. [5]) assigns 
E nr = 8.7 ± 2.7 keV to event A, and E„ r = 6.9 ± 2.8 keV 
to event B, in recognition of the smaller electron signal. 

As shown in Fig. [T] the fraction of energy given to 
electrons (/„) is expected to (and in fact appears to) de- 
crease with decreasing E nr . But the rate of decrease of 
C e tf does not necessarily follow the same slope. This 
is a consequence of recombination, as shown in Fig. [2] 
As of this writing, a direct calibration of the available 
signal from nuclear recoils with energies E nr < 4 keV 



has remained out of reach. Yet broad-spectrum neutron 
calibration data [9] , indicate that liquid xenon direct 
detection experiments have achieved a sensitivity to nu- 
clear recoils of lower energy. The question of "how much 
lower?" has necessarily relied on extrapolations of C e ff- 
In the absence of any theoretical expectation, a reason- 
able choice has been to assume that C e ff — 0.115 (flat) 
for E nr < 4 keV [2J. We have already shown that this is 
incompatible with theoretical expectations, if k = 0.110. 
It remains so even if we optimistically assume that Lind- 
hard's calculation of k = 0.166 is correct. For example, 
the expectation at E Jlr = 2 keV is then C e ff — 0.036 
for \-T B = 0.2, and £ eff = 0.089 for 1 - T e = 0.5. 
Given that in all reasonable cases, theoretical expecta- 
tions suggest that C e ff should decrease with decreasing 
E,„. , it is difficult to support the more aggressive of the 
two exclusion limits (solid curve, Fig. 5 of [2]) proposed 
by the XENON100 collaboration. 
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